LU TP 12-26 
20 June 2012 



The Anatomy of the Pion Loop Hadronic Light by 
Light Scattering Contribution to the Muon Magnetic 

Anomaly 



Mehran Zahiri Abyaneh 

Thesis advisor: Johan Bijnens 

Department of Theoretical Physics, Lund University 
Solvegatan 14A, S22362 Lund, Sweden 

Abstract 

This thesis investigates the Hadronic Light by Light (HLL) scattering contribution 
to the muon g — 2, which is one of the most important low energy hadronic effects 
and consists mainly of the quark loop, the pion pole and the charged pion and 
kaon loops. In this work the charged pion loop has been investigated more closely. 
After reviewing the subject a preliminary introduction to Chiral Perturbation Theory 
(ChPT), Hidden Local Symmetry (HLS) model and the full Vector Meson Dominance 
(VMD) model is given, and they are used to calculate the pion loop HLL scattering 
contribution to the muon anomalous magnetic moment. The momentum regions 
where the contributions of the bare pion loop, the VMD model, and the HLS come 
from, have been studied, to understand why different models give very different 
results. The effects of pion polarizability and charge radius on the HLL scattering, 
which appear at order p 4 in ChPT, from Lg and L\q Lagrangian terms and their 
momentum regions have been studied. 
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1 Introduction 



1.1 Theory 

Elementary particles have some inherent properties including charge, mass, spin and life- 
time. As important as these quantities, are the magnetic and electric dipole moments 
which are typical for charged particles with spin. Classically, an orbiting particle with 
electric charge e carrying mass m entails a magnetic dipole moment given by 

M=£L, (1.1) 

where L is the angular momentum of the particle. Magnetic and electric moments interact 
with external magnetic and electric fields via the Hamiltonian 

H — —fjL • B — d • E , (1.2) 

where B and E are the magnetic and electric field strengths and \x and d the magnetic 
and electric dipole moment operators. The magnetic moment is often measured in units 
of the Bohr magneton /j, B which is defined as 

a B = — = 5.788381804(39) x lO^MeVT" 1 , (1.3) 
where T stands for Tesla. When it comes to spinning particles, the angular momentum 



operator in (1.1) should be replaced by the spin operator [TJ. 

For a charged elementary particle with intrinsic spin and charge q, the magnetic moment 
is written 

M = 9 .^S, (1.4) 

where, S is the spin operator. The constant g s is the Lande g-factor. Although the Dirac 
equation predicts that g s = 2 for electron-like particles, it is slightly greater than 2, and 
theoretically it is useful to break the magnetic moment into two pieces 

fi=(l + a)— , (1.5) 
m 

where a = The first piece, called the Dirac moment, is 2 in units of the Bohr magnetic 
moment. The second piece is called the anomalous (Pauli) moment, and a is a dimensionless 
quantity referred to as the anomaly. 

In 1947, Schwinger, having managed to eliminate divergencies arising in the calculation 
in loop corrections in QED, showed that the deviation of g s from 2 can be ascribed to 
radiative corrections. The first order correction known as the one-loop correction to g = 
2, is shown diagrammatically in Figure [TJ More generally, the Standard- Model corrections 
to the electron, muon or tau anomaly, a(SM), arise from virtual leptons, hadrons, gauge 
bosons and the Higgs boson. This includes the dominant QED terms, which contain only 
leptons and photons; terms which involve hadrons including hadronic vacuum polarization 
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Figure 1: The Feynman graphs for: (a) g = 2; (b) the lowest-order radiative correction 
first calculated by Schwinger. Figure from [2J. 

and hadronic light by light (HLL) corrections, and electroweak terms, which contain the 
Higgs, W and Z. That is, the anomaly for lepton I is calculated as 

ai = a® ED + af ad + aY eak . (1.6) 

An introduction to the theory can be found in [3]. It should be mentioned that, in the 
Standard Model calculations of ai, all contributions coming from the mass scale mi ^> M 
in loops are suppressed by powers of mi/M, and all with in the range M 3> mi are 
enhanced by powers of ln(mi/M). Therefore, for the electron, the most important parts 
come from the QED part where the mediator is the massless photon [1] and the sensitivity 
to hadronic and weak effects as well as the sensitivity to physics beyond the SM is very 
small. Typical Feynman diagrams which contribute to the electron magnetic anomaly are 
shown in Figure [2j This allows for a very precise and model independent prediction of a e 
and hence to determine the fine structure constant a with the highest accuracy, which is 
needed as an input to be able to make precise predictions for other observables like a M . 
This could be done, matching the predicted value of af M [3] 

af M = 0.5--0.32847844400(-) 2 + 1.181234017(-) 3 -1.7502(384)(-) 4 + 1.70(3)x 1(T 12 , 

71 71 71 71 

(1.7) 

where the hadronic and weak contributions are also accounted for, with the observed value 
afP = 0.0011596521883(42) to find [3] 

ar l {a e ) = 137.03599875(52) . (1.8) 

This value is six times more accurate than the other best assessment via the quantum Hall 
effect, which returns 

a~\qH) = 137.03600300(270) . (1.9) 

As discussed above, the QED contributions to are the same as for the electron 
however, the heavy leptons are also allowed inside the loop this time. The overall QED 
contribution to a M then reads jl] 

oed = 11658471.809(0.016) x lO" 10 . (1.10) 
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Figure 2: Typical second and third order QED loop corrections. Figure from [2j. 

On the other hand, is much more sensitive to all three types of effects accounted above, 
and even to physics beyond the Standard Model due to the higher mass of the muon [lj [2]. 

The Electroweak contribution to a M is divided into two parts, one and two-loop contri- 
butions as shown in Figure [3j so that 



,EW 



+ < (2) 



;i.ii) 



which results in 



a 



EW(1) 
M 

EW(2) 
M 



19.48 x KT 10 
-4.07(0.1)(0.18) x KT 10 
15.4(0.1)(0.2) x KT 10 . 



;i.i2) 



Both the QED and electroweak contributions can be calculated to high precision. In 
contrast, the hadronic contribution to cannot be accurately evaluated from low-energy 
quantum chromodynamics (QCD), and leads to the dominant theoretical uncertainty on 
the Standard- Model prediction [2]. In fact, since effects of the energies higher than the 
muon mass are suppressed by powers of (m^/M), the relevant QCD contributions to a M are 
in the non perturbative regime. Nevertheless, there exists a consistent theory to control 
strong interaction dynamics at very low energies, which is called chiral perturbation theory 



(ChPT) [5] and will be discussed in Sec. 2.4 



The hadronic contribution is divided in two parts: the hadronic vacuum polarization 
contribution Figure [4], and the HLL, Figure |5j that is 
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Figure 3: Electroweak one loop and two loop contributions to a M . Figure from [2J. 




Figure 4: The hadronic vacuum polarization contribution, lowest and higher orders. Figure 
from [2] 




Figure 5: Hadronic light by light contribution. Figure from [6]. 
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a had = a (hvp) + a (HLL) _ ( 113 ) 

The vacuum polarization is divided into the leading order and next-to-leading order, whose 
contributions are [7j 

a Had,LO = 6 9o.9(4.4) x KT 10 (1.14) 

and 

a Had, HO = _ 9 8 (o.l) x 1(T 10 . (1.15) 

The part we are interested in in this work, is the hadronic light by light scattering, which, 
contrary to the vacuum polarization part, can not be expressed fully in terms of any 
experimental data and should be dealt with only theoretically and hence, it can be a 
source of more serious errors [TJ and makes the result model dependent. It consists of three 
contributions, the quark loop, the pion exchange and the charged pion (Kaon) loop [2]. 
Due to considerations of the Ref. [8], the estimation of the HLL contribution to the muon 
g — 2 is 

a (h.LxL) = (la5 ± 2 g) x 1Q -10 ? ( L16) 

which is suffering from a large error, as discussed above. 

Calculating the HLL part is the trickiest. Although, ChPT is a reliable theory of 
hadrons at low energies, its usage for the pion exchange brings about divergences and 
one should resort to certain models to get rid of them. One can just introduce some cut 
off energy, but, the way to do it systematically is to cover the photon legs with vector 
mesons. These vector mesons cure the infinities similar to what the Pauli Villars method 
does in QFT, although the Pauli Villars is a pure mathematical manipulation, while vector 
mesons are observable physical entities. There are certain models to do the job (below). 
Historically, after Ref. [9] calculated the HLL part via the naive VMD approach, which 
does not obviously respect the electromagnetic Ward identities [10J, the first thorough 
consideration, compatible with the Ward identities, was by Bijnens, Pallante and Prades 
[TT| IT2"] via the Extended Nambu-Jona-Lasinio approach, assuming full VMD. The other 
was by Hayakawa, Kinoshita and Sanda [T3] using the HLS model. Then, Knecht-Nyffeler 
recalculated the 7Tq,t],t]/ exchange contribution via the quark-hadron duality in the large 
N c limit of QCD [TJ], and found a sign difference with the previous results. Subsequently 
authors of both previous works found a sign mistake which was corrected [15j. Meanwhile, 
afterwards, matching between the short and the long distance behavior of the light-by-light 
scattering amplitude, Melnikov and Vainshtein found some corrections [16J. 

However, as mentioned above, the HLL contribution consists of three parts among 
which, we are interested in the charged pion loop correction in this work. The reason is, as 
can be seen from Table [T], different approaches to this part led to very different results. In 
fact, when the vector mesons are introduced into the calculation, one expects that results 
are heavily suppressed, compared to the bare pion loop case. However, as both VMD and 
HLS models use this mechanism, one might wonder, why the full VMD result is about 
three times larger than the one from the HLS one. This is the main question which is tried 
to be answered in this work. 
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Figure 6: Spin precession in the g — 2 ring. Figure from pQ. 



1.2 Experiment 

A diagrammatic scheme of the a M measurement is shown in Figure [6] pQ. To measure the 
magnetic anomaly an electric field E and/or a magnetic field B must be applied. The 
general formula, derived by Michel and Telegdi [3] in 1959 for this purpose, reads 

e f „ 

u> a = u> s - u> c = <^ a„B - 

where ui c = eB/m^ is the cyclotron frequency, u s = eB /m^rf + a^eB /m M , 7 = l/y/l — v 2 
and v the muon speed. If one forgets about the electric dipole moment of the muon,<i M , 
so that u a is independent of d^, and chooses 7 such that a M — l/(7 2 — 1) = 0, which 
corresponds to the energy 3.1 GeV, called the magical energy, the measurement of a M 
reduces to measuring the magnetic field and the value of u a . As for co a , one should note 
that the direction of the muon spin is determined by detecting the electrons resulting from 
the decay /i - — > e~ + v e + z/ M , or positrons from the decay of /i + as shown in Figure [7j 
The number of electrons detected with an energy above some threshold E t , decreases 
exponentially with time as shown in Figure [8j according to the formula 



Charged pion and Kaon Loop Contributions 


a M x 10 10 


Bijnens, Pallante and Prades(Full VMD) 
Hayakawa and Kinoshita (HGS) 
Kinoshita, Nizic and Okamoto (Naive VMD) 
Kinoshita, Nizic and Okamoto(Scalar QED) 


-1.9 ±0.5 
-0.45 ±0.85 
-1.56 ±0.23 

-5.47 ±4.6 



Table 1: Results of different approaches to the charged pion loop HLL contribution to 
a M U [13] ■ 



a u + 



r 



2d 

(3xE\--^{(3xB + E} , (1.17) 
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Figure 7: Decay of ji + and detection of the emitted e + . Figure from pQ. 
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Figure 8: Distribution of counts versus time. Figure from [T]. 
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N e (t) = N G {E t )e- thT - {1 + A{E t ) cos[(u a )t + <£>(E t )}} , 



(1.18) 



where r M is the muon's lifetime in the laboratory frame. This allows one to extract u a . 
Then, one uses the relation 

B = £-, (1.19) 

between the Larmor spin precession angular velocity of the proton, u p , the proton Bohr 
magneton, fi p , and the magnetic field B, to obtain 



a ^=A^R' (L20) 



where R = u a /u p and A = /i M //i p with //^ the muon Bohr magneton. The value of A is 



measured separately and is used by the experiment to obtain via the relation (1.20). 

Before the E821 experiment at Brookhaven national laboratory between 2001 and 
2004 PP, results of a series of measurements accomplished in the Muon Storage Ring at 
CERN were in good agreement with theoretical predictions of the Standard Model of par- 
ticle physics, that is 

a ex P = 1165924.0(8.5) x 10" 9 ajf = 1165921(8.3) x 10~ 9 . (1.21) 

The BNL experiment managed to improve the CERN experiment 14 fold. The BNL average 
value is [T7] 

a M = 11659208.0(3.3) [6.3] x 10~ 10 , (1.22) 

where the uncertainties are statistical and systematic. The comparison between the exper- 
imental and theoretical values has been done in Figure [9] As can be seen, judging by the 
experimental accuracy achieved in the past decade at BNL, a small discrepancy at the 2 
to 3 a level has persisted with the theoretical predictions. This discrepancy is still debated 
and many conjectures have been made to link it with physics beyond the standard model. 



1.3 Overview of this work 

However, as mentioned above, the theoretical predictions in the realm of the SM are still 
obscured by the hadronic calculations. This work will try to address the charged pion 
loop, as a part of the HLL scattering contribution to a M . The structure of this work is as 
follows. In Sec. |2j QCD and its chiral symmetry will be discussed to give an introduction 
to ChPT. Sec. [3] is devoted to the Hidden Local Symmetry (HLS) model, as an extension 
of the ChPT. Sec. [4] will have a closer look into generalized Feynman vertices for different 
models and some short distance constraints. In Sec. [5] the main body of calculation of 
via different models is discussed and Sec. [6] deals with the role of different momentum 
regions contribution to a M . Finally, in Sec. [7| conclusions and prospects are given. 
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Figure 9: Comparison between the theoretical and experimental values of the a M , experi- 
mental results in top and theoretical values in below. Table from [IS] . 

2 QCD and chiral Symmetry 
2.1 Effective field theory 



There is a folklore theorem ascribed to Weinberg which states [19J: For a given set of 



asymptotic states, perturbation theory with the most general Lagrangian containing all 
the terms allowed by the assumed symmetries will yield the most general S-matrix ele- 
ments consistent with analyticity, perturbative unitarity, cluster decomposition and the 
assumed symmetries. In other words, regardless of the underlying theory, when the de- 
grees of freedom and the symmetries relevant to the energy scale at hand are known, the 
effective Lagrangian built based on them will address the same physics as the underlying 
theory. If a small parameter, A is also realized in the effective theory, one can conduct 
perturbative calculations upon this parameter. Having this in mind, one would go ahead 
with constructing an effective theory of strong interactions in low energies where the orig- 
inal QCD Lagrangian runs into problems due to the fact that in the regime p 2 < 1 GeV 
squared, where the meson dynamics take place, the QCD coupling constant is large. In 
this energy regime, the fundamental particles are hadrons rather than quarks and gluons. 
To build a Lagrangian for a process happening at a scale p C A, one can use a expansion 
in powers of p/A where A is the cut-off energy of the model. Then, the Lagrangian could 
be organized as a series of growing powers of momenta, i.e. of derivatives as 

C = C 2 + C 4 + --- + C 2n + -- - , (2.1) 

where the subscript indicates the number of derivatives. After building the Lagrangian 
like this, one should use Weinberg power counting to realize to which order a given diagram 
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Figure 10: Left is the interaction, mediated via the a particle, right is the same interaction 
when the mediator has been integrated out. 



belongs. Based on the Weinberg power counting scheme, the most important contribution 
to a given scattering comes from the tree level diagram £2- The contribution from one 
loop diagrams is suppressed with respect to tree level and is the same size as the level 
of contribution from the lagrangian £ 4 . The formula determining the power counting 
quantitatively for a given diagram is 

00 

D = 2 + 2(n - l)N 2n + 2N L , (2.2) 

71=1 

where N2 n is the number of vertices originating from £2™ and Nl is the number of loops. 



2.2 Linear sigma model 

To understand how an effective theory can describe dynamics correctly, forgetting about 
the underlying theory, we digress to the linear sigma model as an example. Let's start with 
the linear sigma model Lagrangian 

L = ^0-^0 T -/x(0-0-a 2 ) 2 , (2.3) 

Where the vector field 4> = (0i, • • -,<Pn) is a N-component real scalar field. The field 
has a nonzero vacuum expectation value that is 0q = 4>i + ' ' ' + 4>n = 0?. We assume 
that among an infinite number of ground states that satisfy this condition, one of them 
is chosen dynamically so that, the symmetry is spontaneously broken to the sub group 
H = 0(N — 1). This leads to generation of iV — 1 Goldstone bosons according to the 
Goldstone theorem [20], which are taken to be tt\ Taking a simple choice 4> = (a, 0, ..,0) 
and expanding cf> around O , assuming ti C a as the parameter of expansion and then 
integrating out the a field, one finds the corresponding effective Lagrangian of the non 



linear sigma model. In Figure 10 the difference between the case when the sigma particle 



exists and when it is integrated out in the limit of p /m a <C 1 is depicted. 



2.3 Chiral symmetry 

Returning to QCD, one can observe that the degrees of freedom to be dealt with at low 
energy, namely baryons and mesons are quite different from quarks and gluons which are 
the main players at high energies and the method described for the sigma model to integrate 
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out the heavy field and build up an effective theory does not seem to be applicable here. 
However, the fact that there exists an energy gap between the family of pseudo scalar 
mesons and the rest of the hadrons, and that they can be accounted for as Goldstone 
bosons of a broken symmetry, encourages us to look for an effective field theory to describe 
their interactions [5]. The QCD lagrangian is 

L= Hil^ + gsA^-m^i-^G^G^ , (2.4) 

flavors 

where A 11 is the gluon field and G^ v is the gluon field strength tensor with the definition 

G„ v = d„A v - d v A„ - g s fA^A u , (2.5) 

with the / coefficients the structure constants of the group SU(3) co i our . One can define the 
left handed and right handed fields as 

i>R = \(l + 7s)^ ^ = ^(1-75)^ (2-6) 
and ip = ipL + i>R- The QCD Lagrangian when written in terms of ipL and ipR is 

flavors 

+ '$R(i'fd lt + gsA^^ipm - m^LipiR - mi^RipiL - -G^G^ . (2.7) 
If one drops the mass terms the Lagrangian is invariant under following transformations 
ip iL -»■ exp(-kt L ■ X)ipiL ipiR exp(-ia R ■ \)tp iR , (2.8) 

where X a (a = 1,2, • • -,8) are the SU(3) Gell-Mann matrices in the flavor indices. The 
Lagrangian is said to have an approximate symmetry G = SU(3)l x SU(3)r or chiral 
symmetry. Of course quarks are massive and the chiral symmetry is not realized fully 
in nature however, for three lightest quarks u,d,s, it could be assumed to hold approxi- 
mately. But as this symmetry is not visible in the spectrum of light hadrons [5], it should 
be spontaneously broken in nature due to some spontaneous symmetry breaking (SSB) 
mechanism. This leads to the global symmetry G = SU(3)l x SU(3)r to be reduced to 
the subgroup H = SU(3)v- This being the case, the Goldstone theorem dictates that the 
difference between the original number of generators and the final ones, should have turned 
into Goldstone bosons. In the case at hand the number of Goldstone bosons is 8. As the 
chiral symmetry is also broken explicitly due to the quark masses in the QCD Lagrangian, 
the bosons could be recognized as the pseudo scalar mesons, which have acquired a small 
mass due to this explicit symmetry breaking. 
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2.4 ChPT 



2.4.1 Lowest order 



Now that the ground have been laid, one can go ahead by constructing an effective QCD 
theory at low energies. The most general Lagrangian invariant under Lorentz and chiral 
transformations in the lowest order has the form 121 



C 2 = - j-[tr{D^U ] D^U) + 2B tr{UM^ + MU ] )} 



(2.9) 



where F is the pion decay constant in the limits of the massless pion, B is related to the 
chiral quark condensate and U can be shown to be the 577(3) matrix, written in terms of 
the Goldstone fields as 

[7 = exp(iF o 0) , (2.10) 

where 



V 



V2' 1 Ve 1 ! 
K~ 



K° 2 



K+ \ 



The covariant derivative is 
with right and left external fields reducing to 



I, 



-eQA^ 
—eQA^ , 



(2.11) 



(2.12) 



(2.13) 



for this work with e, the electromagnetic coupling and 

l ( 2 
Q=- -1 
6 \ -1 

Now let's see how can one actually calculate with these tools. To find the amplitude 
for the scattering "f(q,s) — > ir(p) + ir(p'), one has 



and hence 



r p. — — e QA^ 



D U U = 3 V U + ieA u [Q,U] 
D v U j = d u U j + ieA u {Q,U j ] . 



(2.14) 



(2.15) 
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Then, starting from the lowest order Lagrangian the corresponding term is 



f 2 TP 2 

— {DJJ(D»U) ] ) = —(d^Ud^U^) 

F 2 

- ieA ll — {Q[(d^UU ] -U ] d^U)}) 

- A„A^([Q,U][Q,tf}) . 



(2.16) 



Putting in from (2.11) and keeping terms only up to second order of 0, the second term 
reads 

C = -e l -A^Q[d^^]) . (2.17) 



Inserting from (2.11) only for the pion field of <fi one gets 



Q[d*<f>, 



/ 2(^7r+7r- - ir+d^-) 

2(<^7r-7r+-7r-<^7r+) 

\ 



and after performing the trace one easily finds 

C = -eiA^TT+Ti- - 7T + & 



Hence, the Feynman rule for the scattering j(q,e) — > it + (p) + n (p'), using 



and 



d 3 p 1 

d 3 p 1 



6 

Yl ( a P £ W ex P(-^ • x) + a'*e r *{p) exp(ip ■ x) 



(2.18) 



(2.19) 



(2.20) 



exp(ip ■ x) + a p exp(— ip ■ x) j , (2.21) 



(2vr)3 ^/%T V 

as the photon and pion field respectively, reads 

M. = iee ■ (p + p) 

and the vertex is proportional to 

Following the same lines for the scattering 7(5, e) + -y(q',e') — > n + (p) + tt~(j/), 
grangian becomes 

C = e%. 

and the amplitude is 



C = e 2 A„A^TT + n 



which leads to the vertex 



M = 2ie 2 e* ■ e 



2ie 2 g flL , . 



(2.22) 

(2.23) 
the La- 

(2.24) 
(2.25) 
(2.26) 
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2.4.2 L 9 and L w 

The Lagrangian in order p A of ChPT has the form 

£ 4 = L^D^D^U) 2 + L 2 {D^D U U){D^D V U) 

+Lz(D»U ] D^UD v U ] D v U) + L A {D»U ] D^U)( X ] U + X U ] ) 
+L 5 (D^D,U( X 1 U + rfx)) + L & ( X ] U + X U ] ? 
+L 7 ( X j U - X U ] ) 2 + L 8 ( X j U X j U + X U ] xU ] ) 
-iL 9 (F^D"UD u U j + F^D^U j D u U) 

+L w {tfF*UF L n, 

where the field strength tensor reads 

F^ R) = 0„J(O* - 9J(r), - i [J(r) M , J(r)„] , 
with 

X = 2B (s + ip) , (2.29) 

in which s, p, and denote the scalar, pseudo scalar, left and right handed external 
fields, respectively (21]. The terms of interest in this Lagrangian for our purpose are 
those containing L 9 and L w . These term correspond to the pion charge radius and pion 
polarizability. 

One can rewrite the term containing L 9 as 

A) = -i(F^ vR D^UD v U^ + F, uL D^D u U) 
= i(D^ uR UD v U ] + D ll F tlvL tfD v U) 



(2.27) 



(2.28) 



+ i{F iluR UD^D u U ] + F^ L U ] D^D U U) = C\ + C 
The second term can be written as 



^9 



: {F^ R U[D^ D U ]U ] + F^lU^D^ D v ]U) 



(2.30) 



(2.31) 



which using the equalities 

[£>„, D V ]W = -i{F^ L tf - rfF^n) [D„, D U ]U = -t(F, uR U - UF^) , (2.32) 

takes the same form as C\q in the above Lagrangian. Since we are dealing only with 
electromagnetic interaction, F^l = F^ uR and C\ becomes 



C\ = {D^ V [UD V U ] + rfD v U\) , 



(2.33) 



Which upon use of the relation (2.15) reads 

Cl = {D^FjudU^ 



ieA v U[Q,tf] 
U ] dvU + ieA u W[Q,U]]) ■ 



(2.34) 
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Then using the definition (2.10), expanding U and keeping terms up to the order 2 one 
finds 



and 



Ud u U ] 
U ] d v U 

U[Q,lP] 
lP[Q,U] 



-d u UU ] 
-d u U ] U 



Fo 



(2.35) 



[Q, 



Using relations (2.11) and (2.12), the final result is 

C\ = -edj^\2(n + d u n- - ^ d v n + 



- 2ieA, y n + - 



7T 



The second part of the Eq. (2.30) can also be calculated in the same way to give 



c 2 



where F, 



flU 



-eQf^ u is assumed and 



f^r = {d^-dyA^A" -d v A») 



(2.36) 

(2.37) 
(2.38) 

(2.39) 



= d ll A v & i A v - d v A^d v A» - d v A^A v + d u A fM d l/ A fl . 

One can also derive the term corresponding to C\q similarly. The only remaining task is to 
derive the part related to and extract the Feynman rules. For example for the 777T7T 
process 

M — (kk'\d fl A v d' J 'A u Tr + ir~\pp'} , (2.40) 
which via using relation (2.20) and ( 2.21[ ), leads to the following invariant amplitude 

M = ie^(p)e>*{p')p^ - Ap^e^e^ip') . (2.41) 

To get a better understanding of what Lg and L w actually represent, the charge radius 
of the pion is related to its electromagnetic form factor in the low energy region via the 
definition 1221 



V 



F (jp 2 ) 

Comparing this with the pion form factor in the low energy limit one gets 



3a 
m 2 n 



(2.42) 



(2.43) 



Meanwhile, it could be shown that |23j, Lg oc 1/m 2 and hence, Lg is proportional to the 
pion charge radius. Using similar arguments, L\q can be shown to be related to the pion 
polarizability. 
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3 Hidden local symmetry model 



As the Hidden Local Symmetry (HLS) is used extensively in this work, we give a brief intro- 
duction to it in this section. In fact, the HLS considers vector mesons as its gauge bosons, 
achieving mass via eating up the Goldstone bosons appearing as a result of breaking of the 
hidden symmetry, which is added to the chiral symmetry of the ChPT Lagrangian [21]. 
Hence this model is a kind of generalization of perturbation theory. Indeed, ChPT at its 
tree level only covers the threshold energy and even after adding the loop corrections, the 
energy it covers is fully below the chiral symmetry breaking energy around lGeV [22]. As 
the energy grows the p meson should be inevitably considered. That is where the HLS 
takes center stage. There are also some other compatible models discussed in the literature 

[22]. 

As discussed above, symmetry of ChPT is of the type Ggi b a i = SU(Nf) L x SU(Nf) R , 
which in the HLS model is extended to G g i b a i X Hiocai with H = SU (Nf)v- It is interesting 
to mention that the HLS model reduces to ChPT in the low energy region when the vector 
meson is integrated out. In the HLS model the variable U of ChPT, introduced in the 



relation (2.10) is divided into two parts 



U 



These new variables can be parameterized as 

£i jr = exp (ia/F a ) exp (±iii/F n ) with 



71 



n a T a , a = a a T a ] , 



(3.1) 



(3.2) 



where n denotes the Goldstone bosons of the global symmetry and have the same definition 
as before, while a denotes those of the local symmetry. These are the Goldstone bosons 
absorbed by the vector mesons to get massive. Also, F n and F a are the corresponding decay 
constants respectively and T a are the generators of the group. Then, one can introduce 
the covariant derivative including the external fields 

D^i = d& - iv^i + iU^ D^ r = d£ r - iv^ r + i&r^ , (3.3) 

where r and I are the same as in (2.13), with the gauge fields of the Hi oca i defined as^] 



v 



1 a 

9 




73 W + ^) 
















<f>n 



satisfying the field strength 



J [LV 



Then, one can build two independent 1-forms out of the above variables 



2i 



a \\{ x ) 



2i 



(3.4) 



(3.5) 



1 This does include an extra £/(l) global but the ui plays no role in this work. 
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Using these 1-forms, one can build the lowest order lagrangian including £j jr and D^i >r to 
the lowest derivative as 

£ = C A + aC v = F>[(c^(x)) 2 ] + F>[(df (x)) 2 ] , (3.6) 

where 

a = FljFl , (3.7) 

is a constant. 

Finally, adding the kinetic term of the gauge bosons, the Lagrangian, in the unitary 
gauge a = 0, takes the form 



C = C A + aC v + CintiVfj) 

= F>[K(x)) 2 ] + F>[(o|(*)) 2 ] - ^ 

= tr[d^d^} + ag 2 F 2 tr[p^] + 2% ( l -ag ) tr [ff[d^ vr] 



2 

- 2eagF 2 AHr[ Pll Q} + 2ie (l - |) A"tr [Q[d„ tt]] 

4 3 r 

+ ae 2 F 2 A^r[QQ] + _^t r [[^[d^Tr]] + ■ • ■ , (3.8) 

where g is the HLS gauge coupling constant, from this one can easily observe that the 
vector meson has acquired mass equal to ag 2 F 2 via the Higgs mechanism. Also other 
couplings could be expressed as 



1 

ig 

2 



9 P = agF v 

g^n = (l - I) e . (3.9) 
The relevant terms of the above Lagrangian for our purpose in this work, are pU] 

^ y ^ y 

tint = -egpA^pl - igpmcp ^ d m tt _ - ig^A^ + d m tt~ 

+ (1 - a)e 2 A f "A^ + 7i- + 2eg p7TW A fl pln + n- . (3.10) 

It should be mentioned that the crucial property of this Lagrangian, regarding our consid- 
eration of the HLL scattering is that it does not have a p°p°7r + 7r~ term. Corresponding 
diagrams to each of the above terms are depicted in Figure 11 Another property of the 
HLS lagrangian is that for the a = 2 case it reduces to the VMD for the pion single photon 
coupling, still it is different from the full VMD version which includes the ppmr vertex as 
well. 
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(a) 2eg imv A» P tl TT+-K- 



(b) -eg p A»p% 




(d) -ig^A^TT+d^ir 
7r-^7r+) 



(e) -ig P ^P fJ .{^ + d' J "K 

Figure 11: Different HLS Lagrangian terms with the corresponding Feynman diagrams. 
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and 777r + 7r vertices 



Up to this point necessary ingredients of a more through discussion of the HLL contribution 
to a M are introduced. As discussed in the introduction, to cure the infinities one has to 
introduce vector mesons in the calculation of the pion exchange, and this could be done 
via VMD models or the HLS model. At this point one is ready to consider what kind of 



change happens to the point diagrams of the ChPT Lagrangian mentioned in the Sec. 2.4 



when the HLS or VMD models are taken into account. In the naive VMD model, one just 
replaces the photon propagator with the term below 

> — : . = — . " „ . (4.1) 



^9pV "VIIA 



r 



r 



mi 



Q 



q 2 — m 2 



However, this simple model is not compatible with Ward identities [10J. 

To proceed more systematically, one can note that the amplitude corresponding to the 
7 to p to 7T7T, depicted in the right of the Figure [12] is 

iee • (p + p' 



M 



a 
>2 



i)(-i)(m 2 g^ - q p q v 



{q 2 - m 2 ) 



(4.2) 



and hence, one only needs to multiply the vertex of jiiti from (2.23) with 



a m 2 p g^ 



9pp 



2 



QpQp 



mi 



(4.3) 



2 q 2 -m 2 p -l q- 

to find the equivalent vertex including both diagrams shown in Figure [12} Following the 



same lines, the amplitude corresponding to Figure 13, is found by multiplying the •yyiiTr 



vertex (2.26) with 



;i - a)9pp9vv ~ 



a f m 2 p gf,p - g M g^ \ 
2 M q 2 -m 2 ) 



(i 



9pp 



9pp9w t 9pp 



ap 9v 



2 

PvPv 



™ 9vv - quqi 



q 2 — m 2 



mi 



p- 



+ 9v 



a q 9pp - q-pqp 



mi 



q 



(4.4) 





Figure 12: The equivalent vertex of the ^tttt in the HLS model. 
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Figure 13: The equivalent vertex of 77^ in the HLS. 



In the full VMD version, one multiplies the point like 77T7T vertex with 

m? — q 2 



(4.5) 



and the jjtttt vertex of (2.26) with the term 



m 2 p g v0 ~ VvVv rrfyw ~ %% 



i 2 p -p 2 



m 2 — q 2 



(4.6) 



These new vertices are fully gauged and chiral invariant as mentioned in Ref . [12] . One 
can also follow the same procedure to retrieve the desired Feynmen rules for the L 9 and 
L 10 . The extension of the point vertex 77T7T is achieved when multiplied with 



(4.7) 



and the amplitude corresponding to the •yyim vertex, including the p 4 corrections, should 
be multiplied with 



9w9vo + 9ixp,Lq (p 2 g U u ~ PpPu) + 9uuL 9 (q 2 g^ - q^) 
+ (L 9 + L w ) (q ■ pg^guu - g^uPfiqv 



where L 9 and L 10 are constants involved in the Lagrangian (2.27). 



(4.8) 



4.1 High energy limit 

In this section, the high energy limits and the matching with low energy limits are consid- 
ered. 
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(a) (b) 
Figure 14: The 777T7T vertex. 

It can be shown that the 777r7r amplitude for two high energy photons with momenta 
Pi ~ —P2 — P is proportional to 1/P 2 . This is done by using the operator product 
expansion for two vector currents and showing that the matrix element of the leading part 
which is proportional to an axial current vanishes. Hence, when P — > 00 the amplitude 
vanishes. The amplitude corresponding to the diagram 14 is 



M 

For P 1 : 



ie 



, f(P lll -2K ltM )(P 2u -2K 2u ) (P 2u - 2K lv ) (J% - 2K, 



V 



2/x J 



(P 2 -K 2 ) z -ml 
P — > 00 this reads 

M = 



29, 



(Pi ~ K 2 ) 



' P 2 



mi 



2g^} . (4.9) 



(4.10) 



which does not vanish fast enough. For the VMD case the above amplitude should be 



multiplied with (4.6). The resulting amplitude vanishes at order P°. Now let us examine 



the HLS case. Then, the first and second term of (4.9) are multiplied with (4.3) and the 



third is multiplied with (|4.4|). In the high energy limit the leading term is 

2 



M 



9nu 



P^Pu 
p2 



(4.11) 



This is only satisfied for a = 1. However, the case HLS with a = 2 does not uphold this 
condition. Hence, one can infer that something must be wrong with it. 

5 Muon magnetic anomaly from light by light ampli- 
tude 

5.1 General 

The response of a muon carrying momentum p to an external electromagnetic field with 
momentum transferred = p — p' is described by the matrix element 



M 



A p u(p')T p (p,p)u(p) 



(5.1) 
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with 



r»(p',p) = F 1 ( P l)Y 



2m, 



F 2 (p 2 3 )a^p 3u -FM)l*° pu V*> + F A (p 2 3 ) [piY-2m lP p 3 } l5 . (5.2) 



The two first form factors are known as the Dirac and the Pauli form factor, respectively. 
In fact [12], the magnetic moment of the fermion in magnetons is p, = 2(Fi(0) + F 2 (0)) 
and in analogy with the classical limit, described in the introduction, one can define the 
gyromagnetic ratio as g = 2/j, and the anomalous magnetic moment as a = (g — 2)/2 = 
F 2 (0) [12J. The form factor F 3 (p 2 ) can be different from zero provided parity and time 
reversal invariance are broken and for F^p 2 ) to be nonzero, parity invariance should be 
broken. Therefore, both are absent in our survey. Since the task of computation of T p (p',p) 
is very involved especially for higher order corrections, one can project out the form factor 
of interest, F 2 (p 2 ) in our case, and then the general form of the contribution can be shown 
to be [3] 

1 -tr[^ + m)r A/3 (0)(^ + m)[7 A , 7/3 ]] . (5.3) 



a 



light— by— light 



48m 



Defining the four point function H pua \ as 



Pl,P2,P3j 



J d 4 xi J d 4 x 2 J d 4 x 3 exp i (pi ■ x x + p 2 ■ x 2 + p 3 ■ x 3 ) 



x (0 I Tj^j^xJUx^jxixs) I 0) 
and using the Ward identity to rewrite it in the form 

5UP^( Pl ,p 2 ,p 3 ) 



U^ x ( Pl ,p 2 ,p 3 ) = -p w - 
the T A ^(0) for the Figure [5] writes 

d 4 pi f d 4 p 2 



Sp, 



3X 



(5.4) 



(5.5) 



T^(p 3 



1 



(2tt) 4 J (2tt) 4 q 2 p\p\{p\ - m 2 ) (p 2 - m 2 ) 
5W^( Pl ,p 2 ,p 3 y 



5p, 



3A 



(5.6) 



with p 4 = p' — p 2 , pa = p — q. The most formidable task ahead is then to build the 



relevant four point functions and to calculate the integral (5.6). One should note that this 



four-point function can be decomposed by using Lorentz covariance as follows 

U^( Pl ,p 2 ,p 3 ) = U 1 (p 1 ,p 2 ,p 3 )g^g^ + U 2 (p l ,p 2 ,p 3 )gP a g^ 

+ U 3 ( Pl ,p 2 ,p 3 )g^g^ 

+ U^ k ( Pl ,p 2 ,p 3 )g^ppi + n 2 ^( Pl ,p 2 ,p 3 )g^pi 

+ ^\pi,p 2 ,p 3 )g p Pp»p a k + ^ k (p u p 2 ,p 3 )g™py k 

+ TL^ k (p 1 ,p 2 ,p 3 )g^ I f j pt + n^ k (p u p 2 ,p 3 )g^p p pl 

+ U^ m ( Pl ,p 2 ,p 3 )p P p^ k p a m , 



(5.7) 
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Figure 15: The pion exchange HLL contribution to a M . Figure from [6]. 



where i } j } k,m = 1, 2 or 3 and repeated indices are summed. There are in total 13£ 
Il-functions. However, due to Ward identities 



PspU p ^{ Pl , P 2, P 3) = q p n p ™P{ Pl ,p 2 ,p 3 ) = , (5.8) 

all of these functions are not independent and using these identities frequently, the overall 
number could be reduced to 43 independent Ii^ krn {pi,p2,p 3 ) functions, 32 of which con- 
tribute to a M , Ref. [12]. When the functions are found, one should add them up, derivate 



them with respect to p 3 , then set p 3 = and put them into the integral (5.6). 



5.2 Integration 

In this work the main focus is to calculate the contribution of the charged pion loop 
light by light scattering to a M . However, to get familiar with the overall idea behind the 
mathematical approach, it would be illuminating to start with the simpler calculation for 



the pion exchange, shown in Figure 15 



5.2.1 Pion exchange 

The key object that is used for this case is the 717*7* amplitude, which can be calculated 
via p] 



d 4 xexpi (p ■ x) x (0 I T{j^(x 1 )j u (x 2 )} | 7r°(p)) 

= ^a/3P>2^° 77 ( m ^Pl'P2) ) ( 5 - 9 ) 

where, Jvtyy is the form factor function and p\ and p 2 are the photon momenta involved 
in the 7r°77 vertex. Calculating this amplitude for each of the vertices of the Figure 15 



constructing the whole amplitude, taking derivative respect to p 3 , putting p 3 = and 



plugging into the relation (5.3) for three different permutations of Figure 15, one finds 
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LbL 



X 



— e 



d 4 pi f d 4 p 2 



(2tt)* 



(27r) 4 plpl(p 1 +P2) 2 [{p + Pi) -m 2 ^[{p-p 2 ) -m 2 ] 
T 2 (pi,P2;q) 



T«o*^ (pI,pI, g 2 )-7v>* 7 * 7 (p|,p|, 0) 



p\ - m 2 



J r 7r o* rr (q 2 ,plpl)J r 7T o* J * 1 (q 2 , g 2 , 0) 
q 2 — ml 



Ti(pi,P2]q 



(5.10) 



with 



16 16 
Ti{pi,p 2 ;q) = y(p-pi)(p-P2)(pi -P2) - -^{P-P2) 2 p\ 

g 

- o (P - Pi) (Pi • P2)P2 + 8(p • P 2 )P?P2 



3 

16 

y 

16 



16 



(p-p 2 )(pi -p 2 ) + ^-«VPiP 2 



— mj(pi -p 2 ) 2 



(5.11) 



16 16 
T 2 (pi,p 2 ;q) = y(p-pi)(p-p 2 )(Pi -P2) - y(p-Pi) 2 P2 

8 8 

+ jj(p • Pl)(Pl ■ P2)P2 + g(P ■ Pl)(Pl ■ P2)P2 
+ | (P • Pi) (P?P2 ) + \ m V\P\ 

- \rnl{p x -p 2 f , (5.12) 

where q — —(pi + p 2 ) has been used in the limit that p^ vanishes. Also, p is the muon 
momentum. This is an eight dimensional integral to be done. In general three of the 
integrations can be done analytically and one is left with a five dimensional integral con- 
sisting of three angles and two moduli. Then, the angles could be reduced to one, using 
the Gegenbauer polynomials technique pQ. Using this technique, the a^ bL can be averaged 
over the direction of the muon in space such that 

<•••>= -L / dfi(p) . (5.13) 



2tt 
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To do so, one defines (4) = (P + Pi) 2 + m 2 and (5) = (P - P 2 ) 2 + 
find pQ 



m 2 with P 



>2 _ 



= —mi, to 



11 1 ( ZX 

arctan 



< (4) (5) > 

< (P-P^ > 

< <™) ^ > 

< W) > 



m 2 R 12 



= (Pi-P2 
= (Pi-P 2 



(1 — R m2 



1 - zt 

2 



(1 — P ml ) 



5m 



(1 — R m l) 

2m 2 

1 (1 - P m2 ) 



< (5) > 



2m 2 



where 



and 



Pmi — \ / 1 + 



4m 2 



0? 



2; 



P1P2 

= 4^( 1_jRml ^ 1_jRm2 ^ 



Also, if: = cos # and # is the azimuthal angle between the momenta P x and P 2 
The integral (5.10) reduces to a three dimensional integral 

(Q 2 + ml) 



3vr 2 



+ 



where 



h{P„P 2 ,t) = 1/iPMQ 2 ' 



(5.14) 

(5.15) 
(5.16) 



(5.17) 



X(P 1; P 2 , f) (8g 2 (P ■ P 2 ) - 2P|(P 2 4 /m 2 - 2P 



- 2P 2 2 Q 2 (2-P 2 2 /m 2 +2(P 1 -P 2 )/mJ)+4P 1 4 

- 4P 2 Q 2 - 2P 2 P|(4 + P 2 /ml - 2P|/m 2 ) + 2/mJ 

- 2Q 2 (1 + (1 - P m i)(P • P 2 )/m 2 ) + P|(2 - (1 - P ml )P 2 2 /m 2 
+ P|Q 2 (1 - R ml )/ml + P 2 (2 + (1 - P ml ) 2 (P • P 2 )/mJ) 



+ 3P 2 P 2 (l-P ml )/m; 



(5.18) 
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pion exchange contribution 


a M x 10 10 


Bijnens, Pallante and Prades 
Hayakawa and Kinoshita [13] 

Knecht and Nyffeler [T^] 
Melnikov and Vainshtein [16J 


5.6 
5.7 
5.8 
7.65 



Table 2: Results of different calculations for the pion exchange contribution to a M . 



I 2 (Pi,P a ,f) = 1/(P 2 P 2 Q 2 ' 



X{P U P 2 , t) (AQ 2 (P 1 ■ P 2 ) + 2P 2 4 - 2P 2 Q 2 + 2P 4 



2P 2 Q 2 - AP 2 P 2 - A/ml) - 2 Q 2 - 3P lQ 2 (l - Rm2)/(2ml) 
3P 2 Q 2 (1 - R ml )l{2ml) + P|(2 + 3(1 - P m2 )P 2 2 /(2m 2 ! )) 
+ (1 - P m2 ) 2 (P • P 2 )/(2m 2 ) + P 2 (2 + 3(1 - P ml )P 2 /(2m 2 )) 



+ (1 - P ml ) 2 (P a ■ P 2 )/(2m 2 ) - P 2 P 2 (2 - R ml - R m2 )/(2m 



The auxiliary function is 



(5.19) 



X{P 1 ,p 2 ,t) 



1 



arctan 



P1P2X 



zx 



1-zt 



(5.20) 



So, using this technique, one manages to reduce the eight dimensional integral (5.10) to a 
three dimensional one. Instead of t one can also use Q 2 = P 2 + P| + 2P\P 2 as variable. 
Results for the different calculations for the pion exchange are shown in Table 2. These 
correspond to different choices of J r 7r o* 7 * 7 (g 2 



,pi,p 2 2 ) 



5.2.2 Bare pion loop 

What was described above for the pion exchange case we now extend to the charged pion 
loop case which is the topic of this work. This is the contribution of scalar QED which 



is renormalizable. For the pion loop after evaluating the four-point function in (5.4) all 



functions in (5.7) are nonzero. 



The contribution of the hadronic light by light scattering arises at the order of (a/iry 



This includes box diagrams shown in Figure [I6J triangle diagrams in Figure [T7j and bulb 
diagrams, shown in Figure 



18 These diagrams, along with their charge conjugates, add 



up to 21. 



Among the contributions to the four point function denoted in the relation (5.7), the 



H i j krn ^p 1 ^p 2 ^p 3 s j ones originate only from the box diagrams of Figure 16, the U^ k (pi,p 2 ,p 3 ) 



ones can originate both from the box and triangle diagrams and the n j (p 1 ,p 2 ,p 3 ) functions 
come only from the bulb diagrams. One needs to find these functions, take their derivative 
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Figure 16: The box diagram contribution to the HLL. 




Figure 17: The triangle diagram contribution to the a M HLL. 
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ft 
Pa 



P2 

Pa 




Pi 



9 




Pi 
h 



Figure 18: The bulb diagram contribution to the a„ HLL. 



with respect to p 3 , set p 3 = and plug into the relation (5.3). To do so, we have used the 
code FORM. 

Let us first illustrate the procedure with the corresponding four point function for the 
first diagram of the Figure [16} This diagram gives 



1 f d d r 



•4 -4 



(2-K) d (r 2 — m 2 )((r + pi) 2 — m 2 )((r + p\ + P2) 2 — m 2 )((r + p\ + p 2 + p 3 ) 2 — m 2 ) 
x(2r + p 1 +p 2 +p 3 )^2r + p 1 ) u (2r + 2p 1 +p 2 ) a (2r + 2p 1 + 2p 2 +p 3 ) f3 . (5.21) 

Using the Feynman parametrization method 

1 C 1 f l ~ x ri—x-y 1 

— — - — 6 dx dy dz— - — - , 

abed J Jq J [a(l — x — y — z) + ox + cy + az] 4 



(5.22) 



one obtains 



6 f d 4 r f 1 , f x ~ x , f x - x - y , 1 

IW/3 = t / 77; v4 / dx / d V dz 
% J (2vr) 4 J J J 

(2r + pi + p 2 + P3) M (2r + Pi)„(2r + 2pi + p 2 ) a (2r + 2p 1 + 2p 2 + p 3 



r 2 — m 2 



P 1 



(5.23) 



with 



m 2 



m 2 + ({x + y + z)pi + (y + z)p 2 + zp^j - xp\ - y(pi + p 2 ) 2 - z(pi + p 2 + p 3 ) 2 
f = r + (x + y + z)pi + (y + z)p 2 + zp 3 . (5.24) 
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To deal with the integrals we have used the relations 



d d r 



(27T) 



f^f v f a fpf(r 



r.2\ 



d d r 



d d r „ „ 2 

(2^p r ^ /(r } 



(2ii) d d(d + 2) 

Ad 1 

u ' 1 ~2 ^ /~2 

(2ttF/ ^ Jlr 



(5.25) 



and that integrals with odd powers in the numerator vanish. The remaining integrals to 
be done are 



d r 



r 4 



d{d + 2)i J (2vr) rf {r 2 - m 2 ) 4 

1 



(r 2 — m 2 ) 2 



1 



+ 



d(d + 2)i J (27r) d 
2m 2 

+ 7^ 



~ 4 

m 



(r 2 — m 2 ) 3 (r 2 — m 2 ) 4 



1 1 A 5 , L 

- + - - 1-Zn— - 1 + -I 



2416tt 2L £ 6 
1 1 A 



24 16tt 2L £ 



1 - In- 



(5.26) 



1 1 y d d f 



d z 7 (2vr) d (f 2 - m 2 ) 4 
1 



dr d 
di J (2tt) c 



1 1 



(f 2 — m 2 ) 3 



1 1 

+ 



+ 



4 167T 2 V 2m 2 6m 2 



and 



d d r 
(2tt) c 



~ 2 

m 



(f 2 — m 2 ) 4 

1 -1 
16tt 2 12m 2 

1 1 1 



ff 2 



m 2 ) 4 



167T 2 6 m 4 



(5.27) 



(5.28) 



where /i is the subtraction scale of the problem at hand, d = 4 — 2s, 1/e = 
l/(e — /n47r + 7 + 1), 7 is the Euler-Mascheroni constant and e creeps in from the di- 
mensional regularization. The expressions for the integral 



d d r 



(2vr) rf (r 2 



m 2 Y 



(5.29) 



can be found in [5]. Then, one evaluates the derivative dU^ap/dps, which means deriving 
m 2 or occurrences of p% in the numerator, and put ps = to find the relevant function to 
be integrated over. A similar procedure is done for the other box diagrams. 
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For the first triangle diagram shown in Figure 17 the four point function reads 



d A r 



i 3 x i 3 



(27r) 4 (r 2 — m 2 )((r + p 2 ) 2 — m 2 )((r + p 2 + p 3 ) 2 — m 2 ) 
x (2r + p 2 ) a {2r + 2p 2 +p 3 ) p , 



which can be parameterized using 



abc 



1 rl-x 



6 / dx 



JO 



dy 



[a(l — x — + 6x + q/] J 



to give 



with 



n 



fiua/3 



m 2 



1 /" d 4 r i 3 x i 3 



—g^ u (2r + p 2 ) a (2r + 2p 2 + p 3 ) 



i J (27r) 4 r 2 — m 

m 2 + ((x + y)p 2 + ?/p 3 j - xp 2 . - y{P2 + P 3 f 
r + (x + y)p 2 + yp 3 . 



Then using again (5.25) and 



11 f d d r 



y*2 



di J (2n) d (f 2 - m 2 ) 3 



1 1 



4167T 2 

1 1 



~ 2 

m 



4167T 2 



1 1 , . 

- H 1 - Zn 

i 2 /i 

1 , m 2 

- - 1 - In— 



2 2 



one gets 



1 f d d r 



1 1 -1 



i J (27i) d {f 2 - m 2 ) 3 416vr 2 2m 2 



And finally, for the first bulb diagram of Figure 18 one has 



d 4 r 



i 2 x i 2 



x2 g^ 2g af} 



(27r) 4 (r 2 — m 2 )((r + p 2 + P3) 2 — m 2 ) 
1 f d 4 r f 1 . 1 

x ^g^QaP 



i J (2ttY 



dx 



167T 2 



~ 2 

m 



--l-ln- 

£ Li 



(r 2 — m 2 ) 2 



(5.30) 



(5.31) 



(5.32) 



(5.33) 



(5.34) 



(5.35) 
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with 

f = r + x(p 2 + Ps) (5.36) 

and 

fn 2 = m 2 + x(x - l)(p2 + P3) ■ (5.37) 

It is important to notice that the divergent parts of each diagrams, when added up, cancel 
each other and the remaining part which contributes to the a M is finite. As a consequence, 
all Inn 2 dependent terms above disappear miraculously. 

After doing the above integrations, depending on the type of the diagram, we apply the 



Gegenbauer polynomial method to perform five of the integrations in (5.6) similar to the 



steps that led to (5.17). The final formula is rather long and is not presented here. Besides 
the Pi,P2,t or Pi,P2,Q integration one is always left with one, two or three Feynman 
parameters that should also be integrated over. One can always shift the parameters in 
the case of box and triangle diagrams so that, the denominator is independent of one 
of them, and it could be analytically integrated out, reducing the size of the expressions 
considerably. It turns out that the different fh 2 for the different box diagrams can all be 
brought in the same form as well reducing the size of the expressions considerably. The 
final integral to be done is a five or four dimensional integral, which we have dealt with 
using the Monte Carlo routine VEGAS. 



5.2.3 HLS 

When trying the same procedure for the HLS case, we can reuse a lot of the previous 



subsection since the vertices are related by (4.3) and (4.4). One should be careful since 



terms including p.r will also appear in the numerator. For example, for the Figure |5j the 
four point function of the HLS is 



n 



fii/a/3 



d 4 r 



i 4 x i A 



(27r) 4 (r 2 — m 2 )((r + pi) 2 — m 2 )((r + p\ + p^) 2 — m 2 )((r + pi + P2 + P3) 2 — rn 2 ) 



(2r + Pl ) v ( g vU - 



P\9vu ~ VxvVu 



p\ — m 2 



(2r + 2 Pl +p 2 )a [9, 



P29aa ~ P2aP2c 



(2r + 2p 1 + 2p 2 + p 3 )p g P p - 



P39f>p-P3()P3p \ 

p\-m 2 J 



(5.38) 
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The four point function of the first diagram of Figure 17 becomes 



n 



d 4 r 



i 3 x i 3 



(27r) 4 (r 2 — m 2 )((r + P2) 2 — m 2 )((r + P2 + P3) 2 — rn?) 



x 9nv 2 ( g^guu + 9[i(i 



a V\9vv 
2 



PwPn 



pi 



(2r + p 2 ) a \9aa ~ 



P29aa ~ P2aP2a 



Pi 



mi 



fn , o , \ /" Pz9pp-PwPw\ 
(2r + 2 P2 +p 3 )^- p ,_ ml j 



and the four point function corresponding to the first bulb of Figure 18 



IS 



<i 4 r 



i 2 x i 2 



(27r) 4 (r 2 — m 2 )((r + P2 + P3) 2 — fn 2 ) 



x 2g liU 2g aP g^g V u + SW 



PwPn 



mi 



pi 



' 9vv 99/ 



mt 



9* 



ap2 2 5- c 



o ~r 9aa „ 9 

' P2 2 m 2 



pi 



J 



(5.39) 



(5.40) 



To deal with terms like p ■ f one can use the relations (5.26), (5.27) and (5.28), but some 
of the fs couple to p. 

Furthermore, as for the infinities in the HLS approach, they only cancel out after taking 
the derivative with respect to p^ and setting p 3 = 0. The HLS is not a renormalizable 
theoryand thus the result could have been divergent but surprisingly, the contribution to 
a M is finite. 



5.2.4 Full VMD 



As it was described in Sec. 4 in the full VMD, •yim vertex is multiplied with (4.5). However 



how it deals with the need of being chiral and gauge invariant, that the naive VMD model 
does not meet. In fact, it could be shown that covering the photon legs with vector 



mesons when using (4.5), because of the Ward identities in 5.8, always the second term 



in this expression cancels and the result is just like multiplying each photon leg with 
m 2 /(m 2 — q 2 ), and this is fully invariant under above mentioned symmetries. 

Now let us illustrate how the four point function changes in this model. Again us- 
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ing (14.51), (14.51) and (15. 211) one finds 



n 



pua/3 



d A r 



i 4 x z 4 



(27r) 4 (r 2 — m 2 )((r + pi) 2 — m 2 )((r + pi + P2) 2 — m 2 )((r + pi + P2 + P3) 2 — 



x(2r +pi +p 2 +p 3 ), 



(2r+pi)„ 



- PiuPiu 



rn. 



pi 



(2r + 2 Pl +p 2 ) a 



9aam 2 p - P 2a P2c 



pi 



(2r + 2 Pl + 2p 2 + p 3 ) 



m 2 -p| 



For the first triangle diagram of Figure 17 the four-point function writes 

^-pua/3 



d 4 r 



i 3 X 2 3 



(27r) 4 (r 2 — m 2 )((r + p 2 ) 2 — m 2 )((r + p 2 + P3) 2 — m2 ) 



x 



V m 2 — o? m 2 — a 2 / 



(2r + p 2 ) ( 



P Pi 

2 



P2aP2c 



m 2 -p\ 



(2r + 2p 2 +p 3 ){ 



m;-p w p^ \ 

m 2 — vi ) 



and for the first bulb diagram of Figure 18 one has 



i 2 x i 2 



(27r) 4 (r 2 — m 2 )((r + p 2 + p 3 ) 2 — m 2 ) 

_ / mj^* - PipPiu m 2 p9pji - q^qp \ 
x 2g^zg al 3 [ — — ^ ~2 — ~2 — J 



m 2 — p\ 



m 2 — q 2 



x 



m 2 -p\ 



m 2 p -p\ 



1 ■ 



(5.41) 



(5.42) 



(5.43) 



5.2.5 Lq and L 



10 



Since Ref. [25] argued that order p 4 effects might be important, we have also calculated 
the contributions of L$ and Lio to this order. Using the previous results, the four point 
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function of the Figure [5j taking into account the Lg and Lio corrections, takes the form 

11- fiua/S 



d A r 



i 4 x i A 



(27r) 4 (r 2 — m 2 )((r + pi) 2 — m 2 )((r + pi + P2) 2 — m 2 )((r + p± + P2 + P3) 2 — m 2 ) 
x (2r +P1+P2+ p 3 )n (gpji + L 9 (q 2 g^ - q^) 

{2r+p 1 ) v [g vD + L 9 (pjg uP - p lu p w )^j 

(2r + 2pi + p 2 ) a (g a a + L Q (plg a a - P2aP2c 

(2r + 2pi + 2p 2 + P3)p (f)pp + Lg(j%gpp - PsfsPsp) 
The same function for the first triangle diagram of Figure 17 writes 



n 



li.va.fi 



d 4 r 



i 3 x i 3 



(27r) 4 (r 2 — m 2 )((r + P2) 2 — m 2 )((r + £>2 + P3) 2 ~ m 2 ) 
x V (gw9vi> + 9^ L 9 (plgw ~ PwPiu) + 9uuL 9 (q 2 g^ - q^) 

+ (L 9 + L 10 ) (q ■ pig^guu - g^uPifiqc 

(2r + p 2 ) a (g aa + L 9 (plg 

act P2aP2a) J 

(2r + 2p 2 + p 3 )p (g p p + L 9 {p\g^ - pspPsp) 



and the four point function of the first bulb diagram of the Figure 18 



IS 



d 4 r 



i 2 x i 2 



x 



(5.44) 



(5.45) 



(27r) 4 (r 2 — m 2 )((r + p 2 + p^) 2 — m 2 ) 
2g^2g aP (g^g uD + g m L 9 (p\g uD - PwPi u ) + guuL 9 (q 2 g^ - q^) 

+ (L 9 + L 10 ) (q ■ pxg^gvv - g^p^t 

9aa9ft + 9aaL 9 {$9 ft ~ P^Pw) + gtfLg {p^aa ~ P2aP2a) 

+ (Lg + Lio) (p 2 ■ Pzgaagft ~ 9a0PvslPz^) ) ■ (5-46) 

Methods for doing the integrals are the same as of the HLS part. Furthermore, in this 
case the n /iI/Q/3 are finite whereas, contribution to not finite. 
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Cut off 

v lit UU 


io 10 ^ 


GeV 


bare 


VMD 


HLS a = 2 




HLS a = 1 


0.5 


-1.71(7) 


-1.16(3) 


-1.05(1) 


-1.64(1) 


-1.35(1) 


0.6 


-2.03(8) 


-1.41(4) 


-1.15(1) 


-1.80(1) 


-1.59(1) 


0.7 


-2.41(9) 


-1.46(4) 


-1.17(1) 


-1.85(1) 


-1.76(1) 


0.8 


-2.64(9) 


-1.57(6) 


-1.16(1) 


-1.79(1) 


-1.88(1) 


1.0 


-2.97(12) 


-1.59(15) 


-1.07(1) 


—1.53(1) 


-2.03(1) 


2.0 


-3.82(18) 


-1.70(7) 


-0.68(1) 


+ 1.15(1) 


-2.16(1) 


4.0 


-4.12(18) 


-1.66(6) 


-0.50(1) 


+6.18(1) 


-2.14(1) 



Table 3: Results of different calculations for the pion loop. 

6 Relevant Momentum Regions for the pion Loop 
Contribution. 

6.1 Dependence on the photon cut— off A 

Up to this point, the process needed to be taken to calculate the a M for each approach 
has been described. It is now relevant to consider the way a M behaves under change of 
momenta. In fact, one expects that since the HLS model and the ChPT in higher orders 
are non renormalizable, it should be somehow visible through the a M as well. Furthermore, 
after all one of the main goals of the present work has been to deal with differences of the 
VMD and HLS model. One can also expect that the differences should somehow reveal 
themselves via these considerations. To see how, we have calculated different values of 
for different cut-offs so that Pi <A,P 2 < A,Q < A and results are shown in Tabl^3j 

As can be seen, the HLS and the VMD have the same behavior in low momentum 
region however, for higher momenta, HLS starts to behave oddly. Meanwhile, the Lg,Lio 
behavior for higher momenta is as expected due to the non-renormalizability of ChPT in 
order p 4 . 

6.2 Anatomy of the relevant momentum regions for the pion 
Loop Contribution. 

But, as it has been discussed in the Ref. [6], one can also investigate how different regions 
of momentum, contribute differently to the HLL contribution to the a M . The technique 
is the same with calculating the whole value of but this time, instead of taking the 
integral over all variables, one leaves P\ = —pi, P$ = —v\ an d Q 2 = —q 2 unintegrated. 
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Figure 19: of the Eq. (19) as a function of Pi and P 2 for the VMD choice for the pion 
exchange. a M is directly related to the volume under the surface. Figure from [6]. 



Mathematically we write 



dhdl 2 dl q a^ LQ (l u l 2 ,l q ) 



(6.1) 



with li = log(Pi/GeV), l 2 = log (P 2 / GeV) and l q = log(Q/GeV). To see exactly how the 
momentum region above the 1 GeV contribute to the a M , it is better to use logarithmic 
scale as has been discussed in [gj. Also, the total amount of is proportional to the 



volume under the surface of each diagram. For example, as is shown in the Figure 19, the 



most important contribution of the pion exchange to the a^, via the VMD model, is coming 
from the low region of momenta, which is expected since, because of the usage of vector 
meson legs, the large momenta are strongly suppressed. Furthermore, the concentration is 
around the equal values of momenta Pi and P 2 . It should also be noticed that, the whole 
value of the a M is proportional to the volume under the surface, after integrating over the 
whole region of momenta. Also, as it is easier to deal with plots with positive values, — a M 
is drawn in the figures for the pion loop. 

As can be seen from Figure [T9| the momentum is concentrated along the line Pi = P 2 . 
Following the same lines as in [6], we have done the same calculation for the Bare, VMD, 
HLS and Lg, L w approaches to the charged pion loop contribution to the a M and compared 
them. For each case we have shown the distribution of a^ L< ^ versus Pi, Q and Q. Figures 



20 



and 21 belong to the bare pion loop case, where the peak is in the low momentum region 
but, a large part comes from the region above 1 GeV. I Figure [20] we also show the cases 



for Pi 7^ P 2 . It is clear that the parts with Pi significantly different from P 2 contribute 
less. In Figure 21 we show the case Pi = P 2 alone. 
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Figure 20: 



,LL 



as a function of different ratios of Pi and Pi versus Q for the bare pion 



loop choice for. — a M is directly related to the volume under the surface. 



Figures 22 shows the VMD case, which is obviously suppressed respect to the bare pion 



loop case, while the pick still lies in the low momentum region. Figure 23 compares the 
bare and the VMD cases. 

2 respectively, where the 



Figures 24 and [25] belong to the HLS case for a = 1 and a 
second one should reproduce the VMD results. Here is where the surprise comes in and, 
as can be seen, the low momentum region peak follows with a dip at the high momentum 
region. The resulting graphes of HLS and VMD are co-plotted in Figure 26, in terms of 
Pi = P2 and Q. Conclusion to be drawn from this diagram is, in both approaches the main 



body of the contribution comes from the low momentum region and in the VMD case the 
large Q tail is larger. But, the large negative contribution to in the HLS side needs some 
justification. This finding gives a better insight into the nature of the difference which has 
led to such a dramatic variation between the full VMD and the HLS results. It should 
be mentioned that, as has been already noticed in the Ref. [TO], the difference would stem 
from the lack of the pfmit vertex in the HLS lagrangian. The case a = 1 in the HLS, which 
has a better higher energy behavior, makes the dip of the a = 2 case vanish. 

and 28 show results of calculation for the Lg, L w and Lg - 



Figures 
and Figure 



27 



29 



— L w cases 

compares them. It should be mentioned that, as ChPT in the order p A is 
nonrenormalizable, the overall value of the a M in these cases are cutoff dependent. There 
is one specific property of the Lg and L\q terms of the p 4 Lagrangian which we would like 
to emphasis that is, when one sets L\q = —Lg, the Lg part contribution behaves like the 
HLS with a = 2 and the VMD part to order p 4 . This could be seen in the Figures 30 



and 31 and could be justified via relation (4.8) for the •yyim vertex, when compared to 



the relation (4.4). Since, resorting to the fact that Lg oc l/m [23], relation (4.8) plays a 



similar role as the relation (4.4) does in the limit L 



no 
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Figure 21: —a^ Q as a function of Pi = P2 and Q for the bare pion loop choice. — a M is 
directly related to the volume under the surface. 




Figure 22: —or ® as a function of P\ = P2 and Q for the VMD choice. — is directly 
related to the volume under the surface. 
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0.1 



71 loop 




Figure 23: —a^ L ® as a function of Pi = Pi and Q for the VMD and the bare pionloop 
choice. — is directly related to the volume under the surface. 



Jt loop 



LLQ 



-Ze-11 
0.1 




Figure 24: —a^ L ® as a function of Pi = Pi and Q for the HLS, a = 1 choice, 
directly related to the volume under the surface. 



-a M is 



40 



71 loop 



HLS a=2 




Figure 25: —a^ LQ as a function of Pi = Pi and Q for the HLS, a = 2 choice. — is 
directly related to the volume under the surface. 



It loop 



HLS a=2 
VIVID 




Figure 26: —a L J jCi of the Eq. (6.1) as a function of Pi = P2 and Q for the VMD and the 



HLS choices. 
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n loop 




L9L10 



-a^ L ® as a function of Pi 



Figure 27: 

related to the volume under the surface. 



P2 and Q for the L 9 , L w choice. 



-a^ is directly 



LLQ 




P] = p 2 



Figure 28: —a^ 1 ® of the Eq. (6.1) as a function of Pi = P2 and Q for the Lg choice. — c 



is directly related to the volume under the surface. 
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ji loop 




Figure 29: — Q as a function of Pi = P 2 and Q for the L 9 and L 9 , L 10 choice. — a M is 
directly related to the volume under the surface. 



JI loop 




Figure 30: —a^ ® as a function of Pi = P 2 and Q for the VMD and L 9 , L 10 choice. — a M is 
directly related to the volume under the surface. 
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K loop 




Figure 31: — as a function of P\ = P2 and Q for the VMD and L 9 choice. — a M is 
directly related to the volume under the surface. 

7 Conclusions and Prospects 

In this work we have recalculated the previous results for the HLL pion loop contribution 
to the muon magnetic anomaly via the sQED, the VMD model and the HLS model and 
all results are in good agreement with the previous ones. To do so, we have extended the 
Gegenbauer polynomial technique to the pion loop case to calculate the integrals involved. 

We have also added the next to leading order ChPT corrections to the lowest order 
results and have shown that, in the corresponding energy region, results are in good agree- 
ment with predictions of other models namely, the VMD and the HLS with a = 2, as 
expected. 

By investigating the momentum regions that each model predicts to have a part in 
the a^ bL , we have found out why the HLS prediction for the pion loop contribution is so 
different with that of the VMD. 

Also, using the OPE approach to the 77717? amplitude, it has been shown that the 
VMD lives up to the expectation but, HLS with a = 2 does not and hence, the HLS can be 
ruled out as a valid model to consider the pion loop contribution to the muon anomalous 
magnetic moment. The a = 1 case which has a better higher energy behavior has final 
results similar to VMD. 

Acknowledgments 

I would like to thank my supervisor for teaching me nearly all I know about this subject. I 
also would like to thank Stefan Lanz for his endless patience with my questions, my family 
in Iran and my friends Behruz Bozorg, Roland Katz, Pouria Jaberi, Donya Mashallah Poor, 
Reza Jafari Jam and Sajjad Sahbaei for their support. 



44 



References 

[1] F. Jegerlehner and A. Nyffeler, Phys. Rept. 477 (2009) 1 [arXiv:0902.3360| [hep-ph]]. 



[2] J. P. Miller E. de Rafael and B. L. Roberts, Rept. Prog. Phys. 70 (2007) 795 |hep- 
|ph/0703049j . 



[3] M. Knecht, Lect. Notes Phys. 629 (2004) 37 |hep-ph/ 030 7239] 



[4] M. Passera, Phys. Rev. D75 (2007) 013002 |hep-ph/ 06 06174] . 
[5] S. Scherer , Adv. Nucl. Phys. 27 (2003) 277 



[6] J. Bijnens and J. Prades, Mod. Phys. Lett. A 22 (2007) 767 hep-ph /0702170 [HEP- 
PH]]. 

[7] K. Nakamura et al. [Particle Data Group Collaboration], J. Phys. G G 37 (2010) 
075021. 

J. Prades, E. de Rafael and A. Vainshtein, (Advanced series on directions in high 
energy physics. 20) |arXiv:0901.0"306| [hep-ph]]. 

[9] T. Kinoshita, B. Nizic and Y. Okamoto, Phys. Rev. Lett. 54 (1984) 717 



[10] M. H ayakawa, T. Kinoshita and A. I. Sanda, Phys. Rev. D 54 (1996) 3137 |hep- 
|ph/9601310] . 

[11] J. Bijnens, E. Pallante and J. Prades, Phys. Rev. Lett. 75 (1995) 1447 [Erratum-ibid. 
75 (1995) 3781] |hep-ph/950525l] . 

[12] J. Bijnens, E. Pallante and J. Prades, Nucl. Phys. B 474 (1996) 379 lhep-ph/9511388] . 



[13] M. Hayakawa, T. Kinoshita and A. I. Sanda, Phys. Rev. Lett. 75 (1995) 790 |hep- 
lph/9503463]. 



[14] M. Knecht and A. Nyffeler, Phys. Rev. D 65 (2002) 073034 lhep-ph/0312"226] . 

[15] J. Bijnens, E. Pallante and J. Prades, Nucl. Phys. B 626 (2002) 410 |hep-ph/0112"255] ; 
M. Hayakawa and T. Kinoshita, Phys. Rev. D 66 (2002) 019902 

Erratum-ibid. D66 (2002) 019902 

[16] K. Melnikov and A. Vainshtein, Phys. Rev. D 70 (2004) 113006 [hep-ph/ 03 12226] . 

[17] G. W. Bennett et al, (Muon (g-2) Collaboration), Phys. Rev. D 7 3 (2006) 072003 
[hep-ex/0602035] . Phys. Rev. D 70 (2004) 113006 [hep-ph/ 03 12226] . 

[18] F. Jegerlehner Acta. Phys. Polon. B 38 (2007) 3021 [arXiv:0703125v3 [hep-ph]]. 

[19] S. Weinberg, Physica A 96, (1979), 327. 

45 



[20] J. Goldstone, A. Salam and S. Weinberg, Phys. Rev. 127, 965 (1962). 

[21] J. Gasser and H. Leutwyler, Ann. Phys. (N.Y.) 158 (1985) 142, Nucl, Phys, B 250, 
(1985), 465. 

[22] M. Harada, K. Yamawaki, Phys. Rept. 381 (2003) 1 

[23] G. Ecker, J. Gasser, A. Pich and E. de Rafael, Nucl. Phys. B 321 (1989) 311. 

[24] M. Bando, T. Kugo, S. Uehara, K. Yamawaki and T. Yanagida, Phys. Rev. Lett. 54 
(1985) 1215 

[25] K. Engel, H. Patel, M. Ramsey-Musolf [arXiv:1201. 0809y2]. 



46 



